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ABSTRACT 


In this paper a cylindrical shell having a very stiff end plate or 
a flange is considered. It is assumed that near the end the cylinder 
contains an axial flaw which may be modeled as a part-through surface 
crack or a through crack. The primary objective is to study the effect 
of the end constraining on the stress intensity factor which is the 
main fracture mechanics parameter. The applied loads acting on the 
cylinder are assumed to be axt symmetric. Thus the crack problem under 
consideration is symmetric with respect to the plane of the crack and 
consequently only the Mode I stress intensity factors are nonzero. With 
this limitation, the general perturbation problem for a cylinder with 
a built-in end containing an axial crack is considered, Reissner's 
shell theory is used to formulate the problem. The part-through crack 
problem is treated by using a line-spring model. In the case of a 
crack tip terminating at the fixed end it is shown that the integral 
equations of the shell problem has the same generalized Cauchy kernel 
as the corresponding plane stress elasticity problem. Even though the 
problem is formulated for a general surface crack profile and arbitrary 
crack surface tractions, the numerical results are obtained only for a 
semi -el 1 iptic part-through axial crack located at the inside or outside 
surface of the cylinder and for internal pressure acting on the cylinder. 
The stress intensity factors are calculated and presented for a rela- 
tively wide range of dimensionless length parameters of the problem. 


1. Introduction 

In recent past, solutions of crack problems in shells proved to be 
quite useful in studying fatigue and fracture of such important structural 
components as pipes, pressurized containers, and a great variety of 
other thin-walled structural elements (see, for example, [l] for applica- 
tions to the fatigue crack propagation of part-through cracks and to the 
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estimation of net-ligament rupture loads in pipelines). The existing 
solutions which are based on either the classical shallow shell theory 
or a Reissner type transverse shear theory have all been given for 
"infinite" shells in the sense that the crack is assumed to be located 
sufficiently far away from the boundaries and all other sources of 
stress disturbance so that all interaction effects may be neglected. 

The differences between the asymptotic crack tip stress vields given 
by the classical theory and by a higher order theory have now been well- 
documented and will not be discussed in this paper (e.g., [2]), except 
to note that, particularly in the presence of local bending, the trans- 
verse shear theory appears to have clear advantages. The problems of 
a cylindrical shell with an axial crack, that with a circumferential 
crack and a spherical shell with a meridional crack have respectively 
been considered ?n [ 3 ], [4] and [ 5 ] by using Reissner' s shell theory. 

The crack problem of toroidal shells with a positive or negative curva- 
ture ratio, including the effects of material orthotropy has been 
studied in [6]. The problem of an arbitrarily oriented crack in a 
cylindrical shell under general loading conditions is considered in [ 7 ]. 
The solutions given in [3]-[7] are all for a through crack. 

The problem of surface cracks in shells is inherently hree- 
dimensional elasticity problem and appears to be analytically intractable. 
There are, however, some numerical solutions based on the technique of 
finite elements [8], [ 9 ], or boundary integral equations [lO]. Recently, 
there has also been some applications of the line spring inode l developed 
in [ll] for plates to surface crack problems in shells (see [12] for the 
results obtained by the transverse shear and [ 13 ] by the classical 
shell theory). Comparison of the results given in [12] with the finite 
element solution given in [9] shows that one could obtain surprisingly 
good results for the stress intensity factors along the border of a 
part-through crack in a shell by using the line spring model with a 
higher order shell theory and with reasonably accurate compliance func- 
tions for the corresponding plane strain edge crack problem under mem- 
brane and bending loads. 
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The primary objective of this paper la to study the Influence of 
a stiffened end on the stress Intensity factors In a pressurized cylin- 
drical shell containing an axial through or a part-through crack near 
the end. Such problems may arise in pipes and cylindrical containers 
having flanges or end plates the bending and membrane stiffnesses of 
which are very high in comparison with those of the shell itself (e.g., 
heat exchanger tubing near the end plates). Thus, in formulating the 
problem it may be assumed that the end of the shell Is '‘fixed", that 
is all components of the displacement and the rotation vectors are zero. 
The part- through crack problem is solved only for a semi -el 1 ipt ic inter- 
nal or external surface crack. However, the technique is quite general 
and can accommodate any crack profile within the confines of basic 
limitations of the line spring model. As in the corresponding plane 
elasticity problems, in the shell problem too the case of the crack tip 
touching the end stiffener requires special consideration with regard 
to the analysis of the crack tip singularity as well as to the method 
of solution. 

2. The Basic Shell Equations 

Referring to Appendix A for normalized and dimensionless quantities 
and to [3]-[5] for details of derivations, in terms of a stress function 
<j>, displacement component w and the auxiliary functions t|> and ft, the 
basic equilibrium equations of a cylindrical shell may be expressed 
as follows: 


^ - <T> W m0 • 

(1) 

V 4 w + X 2 X|(1 -kV 2 ) 0 =• xMl-KV 2 ) 21 p 

(2) 

kV 2 i p - i|/ - w = 0 , 

(3) 

— Ip ■ V 2 ft - ft - 0 , 

(k) 
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where q(x,y) is the transverse shear lending and \p and Si are related 
to rotations as follows [4]: 

a . + *0-^) M 0 « 11 k(1-v) 3.0 

p x ax 2 ay * y ay ‘ 2 ax * 


The shell and crack dimensions and the notation are described in Fig. 1. 
The normalized membrane, moment and transverse shear resultants are 
defined by 


N . ®4 , N -a, N . . -ISL , 
xx ay 2 * yy ax^ xy 3x3y * 


( 6 ) 


33 33 93 33 

— ( ii + V y ) M B -JL- (v ii + 

"xx hX T ' 3x 3y ' * n yy hX^ ' 3x 3y ; 


M = a (13. + ZZ\ 
M xy hxT ' ay riv ' * 


33„ 33^ 

ax 


V a 3 + 3 V «~+3 

v x 3x *x ’ v y 3y p y 


(7) 

( 8 ) 


To solve the problem, one may first consider a cylindrical shell 
without a crack which is fixed at x 2 * 0 plane and which is subjected 
to the given set of external loads. Since the problem is linear, the 
solution of the cracked shell problem may then be obtained by adding 
to this uncracked shell results a perturbation solution obtained from 
the cracked shell with a fixed end by using the equal and opposite of 
the stress and moment resultants from the first solution as the crack 
surface tractions. Thus, in the main crack problem one may assume 
that the transverse shear load q (such as pressure) is zero and the 
crack surface tractions are statically self-equilibrating. 

3. Solution of the Differential Equations 

By eliminating <j> and by assuming that q = 0, from (1) and (2) 
we f i nd 
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V 4 V 4 w + A| (1 »k7 2 ) 



( 9 ) 


In the present study the primary Interest Is In the pressurized cylinder 
problem. Hence, In formulating the problem It will be assumed that the 
plane of the crack Is a plane of symmetry with respect to loading as 
well as geometry. Therefore, for the shallow shell under consideration 
the solution of (9) may be expressed as 


w(x,y) 



fj(x,a)e ' ay dot 



f 2 (y,g)cosgx dg , 


(10) 


Assuming the solution of the ordinary differential equations resulting 
from (10) and (9) of the form 


fj(x,a) * R(a,m)e mx , f 2 (y,g) * S(g,n)e ny , (ll) 

the characteristic equations giving m and n may be obtained as follows: 


m 8 -Aa 2 m 6 + 6a 4 m 4 -(*ia 2 +KX j )a 4 m 2 + a 4 [a 4 +X^(l+K« 2 )] * 0 , (12) 

n 8 -(Ag 2 -heXpn 6 + (6g 4 +KXfg 2 +x!j)n 4 - ^g 6 n 2 + g 8 - 0 . (13) 

Note that, ordered properly, the roots of (12) and (13) have the 
following property 


Re(mj) < 0 , 
Re(nj) < 0 , 


fflj " “fflj * j = l > • • • » 4 » 

O 4 ) 

n. * "Oj , j * 1,...,A . 

(15) 


If we now observe that for the particular shell under consideration 

(Fig. 1) -°°<y<0, because of symmetry it Is sufficient to consider the 

problem for x>0 only, and the external loads are local and are statically 

self-equilibrating (consequently all field quantities vanish as x-*», 

y+-») , by letting R(a,m. )=R. (a) and S (g,n . )*S . (g) , the functions f, and 

J J J J 
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f 2 vanishing respectively at x-« and y— « may be written as 
4 . . m,x .8 n.y 


h m.x 3 

f, (x,ct) » l R. (a)e J , f«(y,3) 

1 j«l J L 


Q ” \ 7 

2 S.(3)e- J , 

J-5 j 


( 16 ) 


where Rj , (j*l , , . . ,4) and Sj, (j*5»*..,8) are unknown functions. 
Similarly, by expressing 

00 00 

<M*.y) * 2 tT 9 j (*> a ) e Jay dot +| g 2 (y,3)cos3x d3 , (17) 


after obtaining w, from the coupled equations (l) and (2) the functions 
g^ <and g 2 may be obtained as 


gj (x,a) 


X jCt 2 4 R, (a) m.x 
“T7 

X 1 Pj 


2 e J » Pj*m 2 -a 2 , (x>0) , 


X 2 8 n 2 S. (3) n.y 


' 4 1 B " j * { \ m / " • 1 

g 2 (y ,3) " tz - 2 -~~i — e J , q »n?-3 2 , (y<0) . 


* s qj 


j j 


(18) 


(19) 


Expressing now ft in the form 


ft(x,y) * 


2tt 


h j ( x,a)e ^ ay da + — 


h 9 (y,8)sin3y d$ , 


( 20 ) 


and assuming that 


h j (x,a) * A(a,r)e r h 2 (y,3) - B(3,s)e sy , (21) 

from (4) it can be shown that the functions hj and h 2 satisfying the 
conditions at x* 00 and y*-« may be expressed as 

r x 

h,(x,ot) » A| (a)e 1 , r, * -[a 2 + (0<x<«) , (22a, b) 


s y 

h 2 (y,3) - B 2 (3)e 2 , s 2 - +[3 2 + -tI— r]* , (— <y<0) , (23a, b) 


k(\ -v) 
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where Aj(a) * A{a,rp and ( 3 ) - B(0,s 2 ) are unknown functions. 

Finally, one may easily show that the remaining differential equa- 
tion (3) will also be satisfied and the solution will have the proper 
behavior at x*», y»-« If It Is assumed that 


>Mx,y) 


2ir 


00 

0 

I 

-00 


(x,a)e ^dct 


2 

IT 


00 

0 2 (y,0)cosex d3 , 
0 


(24) 


0, (x,a) 


4 R. (a) 

r 

1 KPj-I 


,v 


e 9 (y,e) 


B S (3) 
l — 1 — Y 
5 ^j’ 1 


.V 


(25) 


where R . , m . , p . , S., n. and q. are the same quantities which appear 
in (18) and J (l9). J J 

The formulation given above for x>0, y<0 satisfies the conditions 
at infinity and contains ten unknown functions (Rj , S j , Aj and B 2 ). 
These functions are to be determined from the boundary conditions pre- 
scribed on x«0 and y*0. 


4. The Boundary Conditions 

Referring to Fig. I and assuming symmetry with respect to Xj-0 
plane, the boundary conditions for a cylindrical shell fixed at y»0 
and containing a through crack of length 2a along Xj*0, -d<x 2 <-b may be 
expressed as follows (see Appendix A for the normalized quantities) 


u(x,0) = 0, v(x,0) * 0, w(x,0) - 0, 0 (x,0) ■ 0, 0 (x,0) 

x y 

N xy (°,y) * o, M xy (0,y) * 0, V x (0,y) m 0, (-«°<y<0) , 

N X x( +0 *y) “ F,(y> * -d ] <y<-b | , ^ 

u (0 , y ) » 0 , -»<y<-dj, -bj<y<0 , ] 


0, (0<x<°°) , 
(26) 

(27) 


(28a, b) 
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M xx (+°, y ) * F 2 (y ) , -dj<y<-b. , ] 

3 x (0,y) - 0 , -«><y<-d| , -bj<y<0 J 


(29a,b) 


where d| ■» d/a, bj * b/a and Fj and F 2 are the crack surface tractions 

known from the uncracked shell solution. The displacement component 

w Is given by (10) and (16) and the quantities 6 , B , N „ M , V , 

x y xy xy x 

A and may easily be expressed In terms of the unknown functions 
Rj t Sj , Aj and Bj by substituting from the solution given In the pre- 
vious section Into (5) -(8). These expressions may be found In Appendix 
B. The displacements u and v are determined by using the Hooke's law 
and the strain displacement relations 


e U '2 tu i,j +u j > I +Z ,l u 3J +2 J u 3> l ) - <’-J ’ '> z > 


(30) 


where Zlx^^) gives the equation of the middle surface of the shell. 
Observing that Z ^ ^ jj “ “1/R and referring to Appendix A, It 
can be shown that 



2(l+v) 


3y 


ax 


(3D 


* N 

3y yy 


u(x,y) 


- vN , 


obtain 


_L s 

— 00 

(?,+\>)a 2 -m? 

8 (2+v)n? i 

J s-v^ s - i<6)8e 

o J 


(32) 


m.x- lay 


VU.MIK I 


da 


X l 2 

+ yt xw(x,y) , 


( 33 ) 
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-00 J 


7F 


2 f 8 n, n.y 

2 jjy (0 2 +vn|)Sj (0)e J cos0x d$, (x>0, y<0). (3*0 

o J 


By substituting from (33), (3*0, (10), (16b (B.9) and (B, 10) into 
(26) and by inverting the sine and cosine transforms we find 


8 (2+v)n. 2 -$ 2 


M J M j 

2 — r-7 1 S (0) - - 57 

5 9j j 2n 


t 4 ( 2 +\>)a 2 -m . 2 


? p/ r nijW) R j (a)m j d0 


(35) 


8 n.(B 2 +vn 2 ) ,, r°° 4 m 2 +va 2 

sV ' $ J (0) "27 I f R J (a)am J da 

-00 j J 


(36) 


2 S. ( 3 ) 
5 J 


, r* 4 m.R * (a) 

2? I l m7+3 2 ' da » 

'* * i 


(37) 


BS.(B) 


8 p *i vp ' K (i-v) 

5 wij '1 ' ~ 


, f 4 $m.R.(a) 

■2 B 2 (S) " 2iT J f Tkp' - d “ 

00 J J 


K (l-\i) I r «6A, (a) 


r at 

j ~"2la2 " da , 


"Sir j r,*+0 


(38) 


8 n.S (3) I f* 4 am,R. (a) 

l (k Mm) da 

-oo J J 


lU t 

5 < 9 j-l 


k(I-v) f” r l 2A l <a) 


. k(I-v) f r l 

+ TF— j r: 


7 +3 r 


doc . 


( 39 ) 


Consideration of the mixed boundary conditions (28) and ( 29 ) in 
mind, we now define the following new unknown functions: 
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g,(v) 


0 


u(+0-,y) , G 2 (y) - 9y ^ 


pr/ “ 3 y , « 2 \r/ - aU’M* 0 */) » (-«<y<0) . 

From (33), (B,9) and (^0) tt can be shown that 


j F 4 0 +v)« 2 -p, 

(y) * ‘ 5? 3F | f — FT” 1 R J y 

-00 J 


dot 


My) 


j j- 4 nijRj (a) „, ay ^ k(1»v) 


27T 


f 

J 1 KPl- 

-CO J 


00 


(40 


(41) 


f~ ae da " Ti" 


ct^A, (a)e~ lctY da, (42) 


From (28), (29) and (40), If we observe that Gj“0, G 2 «0 for -«<y<-dj, 
-b | <y <C , by defining 


-b, 


Cj (ot) * 


Gj (y>e ,oty dy , (j » 1,2) 


-d 


(43) 


and by using the expressions (B.3), (B,6) and ( 9 , 7 ) it can be shown 
that the homogeneous boundary conditions (27), (4l) and (42) are equiv- 
alent to 


4 m.R. (a) 


it • 

JL_i 
1 Pi 


2 — " 0 • 


4 m • R , (ot ) . 

I " V.-t " fCK(1 “ v)a + a ] C 2 (a) ’ 


(44) 


(45) 


E m.Rj («) - - 1 c 2 (a) , 


(46) 


4 m.R. (a) 2 f 

T 1? j c i (a ' - 


S 

1 Pj 


(47) 


A| (a) - 2C 2 (a) . 


(48) 
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From (44)-{47) one may easily solve for Rj and obtain expressions 
of the form 


Rj (a) - f [Qj (ot)C j (a) + N j («)C 2 (a)] 



2 

E Cj k (a,t)G k (t)dt , (j-I,...,4) 
k® | 


(49) 


Substituting now from (48) and (49) Into (35) - (39) , the remaining 
unknowns may also be expressed as follows: 

-b 

, t r 1 2 

S j (3) - I 2 b k (0,t)G k (t)dt , (j-1 4} , 

■ k« 1 



b 5l ,(M)G k (t)dt 


f 


(50) 


(51) 


where c j<k and bj k are known functions. Thus, once Gj and G 2 are deter- 
mined (48)-(51) and Appendix B would give all the field quantities 
needed. 


5. The Integral Equations 


From the derivations given In the previous section It Is seen that, 
In addition to the assumption Gj-O-Gg for -®<y<-d j ,-b<y<0 used In (43) 

If Gj and G 2 are also assumed to satisfy the conditions 


f l 

J Gj (y)dy ■ 0 , (j-l ,2) , 

-d 

1 


(52) 


Then all boundary conditions (26)-(29) except (28a) and (29a) would 
be satisfied. These two conditions may now be used to determine the 
unknown functions Gj and Gj. Thus, by substituting from equations (B.l), 
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(B.4) and (A8) - (5 1 ) Into (28a) and (29a), we obtain 
-b, 


j » » 

llm E [f V. . (XfCtJe*^ y ^ a da + [ Y. , (y ,t ,3)cosx8d3]G. (t)dt 

x-H-0 1*1 J KJ J KJ - 


- d| J" 1 ~ 


l r k (y)» (k*1 ,2) , (-dj<y<-b j ) 


(53) 


where V^j and Y^j are very complicated but known functions. As t-ty and 
for b*0 as (t,y)-K) the kernels In the Integral equations (53) are expec- 
ted to be unbounded. Since the singular behavior of the unknowns Gj 
and G^ will be dependent on the singular nature of the kernels, It Is 
necessary to examine the asymptotic behavior of these kernels. In (53) 
the Integrands of the Inner Integrals are bounded for all values of 
their arguments. Thus, any singularities the kernels may have must be 
due to the behavior of the Integrands at a * + » and 3 " and these 
singularities may be separated by Isolating the asymptotic parts of 
V^j for |ot|-*» and of Y^j for This can be done In closed form by 

first extracting asymptotic values of m.(a), n.(3), r^ (a) and s 2 (8) 
from the ch^ract&ri'&cfc equations (12), (13), (22b ) and (23b) and then 
by substituting these values Into the expressions of and Y^. . From 
the characteristic equations It can be shown that for large values of 
1 a | and 3 we have 

p. p. 2 

mj (a) * - 1 a | (1 + ~ ^ - g-jr + • ••) » (j**l > • • • >4) , (54) 

n n 2 

(55) 


nj (3) * 3(1 + 2^7 " + • • • ) » (j"5» • • • ,8) , 


r l (a) * “M (1 + kTFvT^ “ 


(56) 


s 2 (3) - 8(1 + ” •••) 


(57) 



« 

-•» 

1 

t 


+ 

n 


f 
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By using (5*0-(57) and by adding and subtracting the asymptotic 
values of V^j and Y kj . In (53), after evaluating the singular terms of 
the kernels In closed form, the Integral equations (53) may finally 
be expressed as follows! 


*b 


f r_L + ^ + *^3 . iL'nl _jl__ + ... yf— ] G ( t ) dt 

j L t-y (v+l)(3-v) t+y 3-v (t+y) 2 3-v (t+y) 2j G l' t; 1 


I 

-d 


r b l 


‘1 


kj , (y ,t)G 1 (t)dt + J k 12 (y ,t)G 2 (t)dt ■ 2 ttF 1 (y) , 

“ d l 

-d 1 <y<-b 1 , 


(58) 


-b. 

/, r 1 . v 2 +6v-3 1 6(1+v) y . Ml+v) y 2 t ^ 

0 v } J C t=y + T^TTTF37 t+y 1=r~ Tt+yF + Tt+yF ] G 2 (t)dt 


»b. -b. 

1 f 1 h 

k 21 (y*t)G^ (t)dt + I k 22 (y,t)G 2 (t)dt - 2^ ~ F 2 (y) , 
-d, -d. 


“ c *l<y < "b| , 


(59) 


where the kernels k j j , (i,j * 1,2) are known bounded functions. It 
should be noted that the Integral equations (58) and (59) must be solved 
under conditions (52). 

For bj>0 the integral equations (58) and (59) have simple Cauchy- 
type kernels and have a solution of the form [14] 

G k <y) " H k (y)/[“(y +d i> (y+hj)]* , -d 1 <y<-b l , (k-1,2) (60) 

where Hj and H 2 are unknown bounded functions. In this case the integral 
equations may be solved numerically by using the technique described, 


" 13 " 


for example, In [ 1 5 ] . On the other hand for bj * 0 (that fs, If the 
crack tip y ■ -b, Is extended to the fixed end), ft Is seen that the 
kernels have singular termlf In addition to (t-y) which become unbounded 
as t+0 and y-K) together. The singular parts of the kernels shown separ- 
ately In (58) and (59) are typically the generalized Cauchy kernels. 
Because of these kernels the solution no longer has the square-root 
singularity given by (60). From (58) and (59) we first observe that 
the dominant kernels in the two equations are identical. Consequently, 
the singular behavior of Gj and G 2 (or u and 3 x ) at the crack tips would 
be identical. This, of course, is the physically expected result. Next, 
to obtain the singular behavior of the solution we express Gj and G 2 as 

G^fy) ■ P k (y)/C-y Y (y+d 1 > aJ ] , 0<Re(y,a))<l , (k- 1,2), (-dj<y<0), (61) 

where the functions Pj and P 2 are bounded and unknown and the unknown con- 
stants y and w are determined by substituting from (61) into (58) and 
(59) and by using the function-theoretic method (see, for example, 

[15]). Thus, from (58), (59) and (6l) the characteristic equations giv- 
ing Y and w may be obtained as follows: 

C0S7TY + ■(■v+lf(3-v) + 2 3-v" V ^ y (y+2 > * 0 * (62) 

COtTTO) ■ 0 . (63) 

The acceptable root of ( 63 ) is 1/2 which is the expected power of 
singularity for the crack tip y*-dj embedded in a homogeneous medium. 

It may be shown that (62) is identical to the plane stress version of 
the characteristic equation for the corresponding plane elasticity prob- 
lem [16]. For example, for v»0,3 it is found that y * 0 .24 165 is the 
acceptable root. In this special case too the system of integral 
equations (58) and (59) may again be solved numerically by following the 
technique described in [153- 


6. Stress Intensity Factors 


After solving the Integral equations the asymptotic behavior of 
the stress state around the crack tips may be examined by using the 
expressions given In Appendix B. For the embedded crack (i.e., for b> 0 ) 
the problem Is the same as in an infinite shell which has been discussed 
previously in detail elsewhere (see, [ 3 ] ~ [6] ) . For the crack tip at the 
fixed end b»0 the asymptotic behavior of the stress state would be 
Identical to the plane stress problem considered in [ l 6] and only the 
related stress intensity factor needs to be given. One should note that 
as shown in [3]"[6], unlike the asymptotic results obtained from the 
classical theory, the membrane and bending resultants given by the 
Reissner's theory have identical behavior around the crack tips and 
this behavior is in turn identical to that given by the plane elasticity 
theory. Thus, to describe the stress state around the crack tips all 
one needs to do is to determine a set of thickness -dependent stress 
intensity factors. 

The in-plane components of the stress state in the shell are given 


by 


a ij ^ X 1 » x 2’ x 3^ = a j j a jj » a |j ( x l» x 2^ 


N. J (x | ,x 2 ) 

, 


& 1 j * *2 * *3 ^ 


12x 

Mjj (x, ,x 2 ) , (f,j»!,2). 




where superscripts m and b refer to membrane and bending stresses, in 
the symmetric problem under consideration only the Mode I stress state 
exists around the crack tips. For the embedded crack the corresponding 
stress intensity factors at the crack tips x 2 a -d and x 2 » -b (Fig. 1 ) 
may, therefore, be defined as follows: 


k.(-b,x-) * lim v / 2(x 2 +bT j (0,x 2 ,x.J , 
x 2+~ b 


k.(-d,xj * lim /-2(x,+d) , (0,x« ,x,) . 

' 5 x 2 +-d L " c i 


(63a, b) 
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These stress Intensity factors may be calculated from the asymptotic 
analysis of the stresses around the crack tips. However, they may 
also be calculated directly from the Integral equations (58) and (59) 
without lengthy asymptotic analysis. To do this we first observe that 
the left hand sides of (58) and (59) give the expressions for N xx 
and M xx on x*0 outside as well as Inside the crack (~dj<y<*bj), and 
only the dominant kernels contribute to the stress singularity. Thus 
referring to (64) and (65), after some calculations It may easily be 
shown that 

kj(-b,x 3 ) * - j ft [H| (-bj) + H 2 (-b ^ ) ] , 

F *7 

kj(-d,Xj) * j /a [H| (-dj) + H 2 (-dj)] , (66a, b) 

where H| and H2 are related to Gj and G2 through (60) and are the 
main calculated results. 

For b « 0 the angular distribution of stresses and the stress 
Intensity factors will depend on the dominant parts of the integral 
equations (58) and (59)* Since the dominant kernels obtained for the 
shell problem are Identical to the plane elasticity problem, (aside 
from the magnitude of the stress Intensity factors) the analysis and 
the results would also be the same. The details of the asymptotic 
analysis may be found in [16]. In this problem It is particularly 
important to note that the stress Intensity factor kj is a measure of 
the amplitude of all in-plane stresses a^., (i,j=l,2) around the crack 
tip y»0, there is no "cleavage" stress in the "second medium" 

(which is rigid) in terms of which kj may be defined, stress state 
around the crack tip may be calculated in terms of kj and the angular 
distributions given in [16], and kj may easily be defined in terms of 
the crack surface displacement (au+x^0 x ) . Similar to [16], from 
the dominant part of the integral equations the stress intensity factor 
at the crack tip may be obtained as 
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kj(0,x 3 ) 



(67) 


After the determination of the thickness-dependent stress Intensity fac- 
tor k| the stress state near the crack tip x 2 »0 may be obtained from 

* ki (o 

a U (r>9>x 3 ) ’ (68) 

where (r,0) are the polar coordinates In (x^.x,) plane and the functions 
f.j may be found in [l6] v< . In this case at the embedded crack tip 
x 2 *-d the relation between the stress Intensity factor and the crack 
surface displacement leading to (66b), namely 

kj(-d,x 3 ) - |> 11m ^ /2(x 2 +d) [u 1 (+0,x 2 ) + x ^, } (+0,x 2 ) ] ( 69 ) 

is still valid and, therefore, referring to (61) kj may be calculated 
from 

X 

k i < ‘ d,x 3 ) '^ [p i ( ‘ d i > + ‘r |, 2 ( ' d i )] • (70) 


7. The Part-Through Crack Problem 

In this paper the part-through crack problem described in Fig. 1 
will be treated by using the line-spring model. The particular version 
of the model used with the Reissner's plate or shell theory is discussed 
in [17] and [12] in detail. Since for the problem under consideration 
the method to be followed and the compliance coefficients to be used 
in the line spring model would be identical to those given in [17] and 
[12], only the results and no analytical details will be discussed in 
this paper. 


FT 


The distributions given for the greatest stiffness ratio in [16] may 
be used as an approximation to the rigid end problem. 
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In the symmetrically loaded shell having a symmetrica! ly oriented 
crack shown In Fig. 1, It Is clear that only the Mode I stress Inten- 
sity factor would be nonzero along the crack front. At a location y 
along the crack front this stress Intensity factor Is calculated from 
[ 12 ] 


k,(y) * /h [ct ( y )g t (g) + m(y)g b (5)] , 5 ” L/H . 


The functions cr and m are defined by 


tf(y) 


N(x 2 ) 

"nr" 



6M(x 2 J 6M(ay) 
h 7 *— » 


(71) 


(72) 


where N and M are the membrane and bending resultants which replace 
the net. ligament under the surface crack and L Is the crack depth 
(Fig. 1). The functions and g b are essentially the shape factors 
for the plane strain problem of a strip of thickness h which contains 
an edge crack of depth L and is subjected to uniform tension or cylin- 
drical bending. These shape factors are given by 

g t (?) - /? (1.1216 + 6. 5200? 2 - 12.38775 4 + 89.0554? 6 

- 188.6080? 8 + 207.3870? 10 - 32.0524? 12 ) , 

g b (S) * St (1.1202 - 1.8872? + 18.0143S 2 - 87.3851? 3 

+ 241.912 4? 4 - 319. 9^02? 5 + 168.0105? 6 ) . (73a, b) 

To determine a and m the integral equations (58) and (59) are 
modified to incorporate the effect of the net ligament. This is done 
by replacing F. and F- (which are equal and opposite to N and M 
calculated from the uncracked shell solution) as follows: 

f, ( y) -*• F,(y) + , f 2 (y) * F 2 (y) + . (74) 
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By using also the linear dependence between the pairs (cr,m) and (u,0 x ) 


[12], the Integral equations are then modified as 

y y -b, 

f - I 1 I 

-Y tt (y) G. (t)dt + Y tb (y) G 2 ^) dt + k (y,t)G, (t)dt 

r . „ 


-d, 


-d, 


" b | 2 

+ 2? I I k Jj (y*t)Gj(t)dt - F|(y), -d 1 <y<-b 1 , 
-d. 


(75) 


+r 


bt 


I 

G,(t)dt - Y bb (y)[ 0 2 (t)dt 


-d 


-d 


a 

lirhl^ 


(y)| 

-d 


-b. 
f I 2 

j J k 2j . (y ,t)Gj (t)dt - F 2 (y ) , -dj<y<-bj , 
-d. 


k s (y,t)G 2 (t)dt 


(76) 


where k g Is the generalized Cauchy kernel given In (58) and (59) » the 
upper (I.e., -) the lower (i.e., +) signs are to be used for the outer 
and the Inner surface crack, respectively. 

After solving (75) and (76) for Gj and G 2 , u and 0 , and then a 
and m are evaluated as follows: 


u(+0,y) 


r y 

G j (t)dt , 3 x (+0,y) 


,y 


d 


1 


G 2 (t)dt 


(77) 


- E Cr tt u + ir tb B x ], m(y) . 6E[y bt u + Y bb e x ] , 


(78) 


where the + and - signs are to be used for the outer and the inner 
crack, respectively. The functions y.j(y), (i,j»t,b) are known in terms 
of g t and g fa (see [12] or [ 17 ]) • 
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8, The Results for a Pressurized Cyljnder 

As an example we consider a cylinder which is fixed in one end 
and is internally pressurized (Fig. 1). The cylinder is assumed to be 
sufficiently long so that the perturbation field of the crack interacts 
with one end only. The crack problem is solved by using the equal ana 
opposite of the stress resultants N 1 j and Mjj obtained in Appendix C 
from the uncracked shell solution as the crack surface tractions in 
the integral equations (58) and (59) or (75) and (76). In the through 
crack problem, as seen from (66), the stress intensity factors are 
linearly dependent on the thickness coordinate Thus, one may dis- 
tinguish a “membrane" and a “bending" component of the stress inten- 
sity factor. We could therefore define the following normalized 
stress intensity factors: 


k (a. ) 
m i 


ki(a,»0) k. (a. ,h/2)-k. (a. ,0) 

(p o R./h)/a b ! (p Q R./h)/a 


(79a, b) 


where oij * -b, identify the crack tip, kj is defined by (65) 

and is calculated from (66), p Q is the internal pressure, 2a is the 
crack length, and Rj is the inner radius and h the thickness of the 
cylinder. The normalizing stress intensity factor (p Q Rj/h ) /a corresponds 
to the value in a flat plate under a uniform membrane stress (p Q R./h) 
arid haying the same size crack as the shell. 

Tables 1 and 2 show the calculated results obtained for certain 
values of the dimensionless length parameters R/h, a/h and c/a of the 
problem, where c defines the crack location. The tables clearly show 
the influence of the end stiffener in reducing the stress intensity 
factors. As c increases the stress intensity factors at both ends 
approach the infinite cylinder values. The effect of the end stiffener 
and the curvature on the stress intensity factors can be seen somewhat 
better in Fig. 2 where some limited information is displayed. The 
figure shows the membrane component of the normalized stress intensity 
factor defined by (79a) as a function of the crack location c for a 
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fixed crack length a"3h and for various curvature ratios R/h. It Is 
Important again to observe that as c-+a, or as the crack tip approaches 
the end stiffener, generally the stress Intensity factors decrease and, 
furthermore, k at x*«-b tends to zero. This may easily be explained 
by noting that for c»a there Is a change In the nature of stress singu- 
larity, l.e., 

k, (-b) 

a a — f or c>a or b>0 , (80) 

U /2r 

k.(0) k.(0)r*' Y 

(7 a for c*a or b*0 , (81) 

,J /i r Y Jlr 

where r is a small distance from the crack tip and kj(-b) and kj (0) are 
finite constants. Thus, from (80) and (81) it Is seen that the stress 
intensity factor based on the definition of the conventional square 
root singularity becomes 

lim k.(-b) - Urn /2r or. . * 1 im k.(0)r^ ^ ■ 0 . (82) 

b+ 0 1 r+0 r-*> 

b+0 

A second somewhat curious observation one may make from Fig. 2 

is that for smaller values of R/h there seems to be a slight "overshoot" 

in the stress intensity factor k (-d) before it tends to the infinite 

m 

cylinder value. The explanation for this may be found in the uncracked 
cylinder solution which indicates that in the perturbation region the 
end stiffener may cause a slight bulging of the cylinder in radial 
di recti on. 

In Tables 1 and 2 the column c/a = 1 corresponds to the crack ter- 
minating at the end stiffener, (i.e., b*»0) . In this case the stress 
intensity factor kj at Xj*!) Is calculated from (67). Despite relatively 
large magnitudes of these stress Intensity factors, as indicated above, 
since the stress singularity at this point is a weaker singularity 


- 21 - 


ORIGINAL PAGE SS 
OF POOR QUALITY 


(l.e., y<£, see (68)), there Is actually a reduction In the amplitude 
of the stress field as b-K). 

The results obtained from the part-through crack solution are given 
In Tables 3“ 1 0 . In all these calculations It Is assumed that the crack 
has a seml-el 1 iptlc profile as shown In Fig. 1 for which the crack 
depth Is given by 


L(x 2 ) * L q /R^F 7? > (-d<x 2 <-b) 


(83) 


or 

L(y) - ^/Hy+e ,) 2 » y ■ x 2 /a, Cj - c/a, (-C| -l<y<-Cj+l ) (8k) 


The stress Intensity factors shown In the Tables are calculated from 
(71) and are normalized as follows: 


k t (x) 


kj (x 2 ) 


ot 




^ g t ^ 0 ) * 5 


o 



kfcW 


k 1 (x 2 ) 


ob 




s b (5 0 ) 



(85) 

( 86 ) 


where x * y+Cj * (x^+cj/a and k fc is the contribution of the membrane com- 
ponent Njj and k fa that of the bending component Mjj of the external 
loading (see, equations (C9 ) and (CIO)). The total Mode I stress inten- 
sity factor aloiig the crack front is then given by 


kj(x) ■ k t (x) + k fa (x) , (- 1 <x< 1 ) . (87) 

The normalizing stress intensity factors k^ and k Qb are the correspond- 
ing edge crack values obtained from the plane strain solution of an 
edge-notched strip under a membrane stress No, - p Q R ,/h or a bending 
moment M<» » p Q R j h/6 . The expressions of g fc and g^ are given by (73). 

For values of L Q /h which are used in the numerical calculations they 
are also given by 


22 


K/ h * s o 

0.2 

0.4 

0,6 

0.8 


1,367^ 

2.1119 

4.035 

11.988 

j 

1,0554 

1.2610 

1.915 

4.691 


In the tables contributions of the membrane and bending loads 
are shown separately in order to give some idea about the nature of the 
loading. The results show that even though generally the membrane com- 
ponent is the dominant stress intensity factor, near the fixed end, 
particularly for relatively small cracks the stress Intensity factor 
resulting from the bending component of the external load may be the 
more significant one. 

Tables 3-7 show the normalized stress intensity factors k t and 
k^ at the maximum penetration point (x * y+Cj - (x 2 +c)/a ■ 0 ) of a 
semielliptlc inner or outer surface crack in the pressurized cylinder 
for various combinations of the dimensionless length parameters l, Q /h, 
a/h, c/a and R/h. In all calculations given in this paper it is assumed 
that the Poisson's ratio v is 0.3. The effect of v on the stress 
intensity factors, however, is known to be rather insignificant [ 4 ]. 

Tables 8-10 show some examples for the distribution of the normal- 
ized stress Intensity factors k £ and k b along the crack front. As 
expected, for cracks very near the end the stress intensity distribu- 
tion is highly nonsymmetric (with respect to the mid-point x * 

(xj+cJ/a * 0 ) and as c increases it becomes mors symmetric (see, for 
example, Table 9, c/a - 10). 

Some sample results showing the variation of the total stress 
intensity factor are given in Figures 3 and 4. The total stress 
intensity factor kj is calculated from ( 85 ) -( 87 ) as follows: 


k,(x 2 ) - k ot k t (J) + k ob k b (J) 


p R. 
-o 1 


HT" g t <V h)+k b^ ),/h/L o 9 b (L o /h)] • (88) 
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Figure 3 shows the stress Intensity factor at the deepest penetra- 
tion point of the crack (i,e,» for L«t Q at X 2 *-c) as a function of the 
distance c to the end of the cylinder for a fixed crack size (a/h ■ I, 

L /h ■ 0,4) and for various values of the curvature. This figure too 
o 

shows the drastic reduction In the stress Intensity factor caused by 
the end stiffener and by the Increase In the cylinder radius. The 
results shown In Fig. 3 are obtained for an outer crack. For greater 
values of R/h If the crack Is very near the cylinder end then the stress 
Intensity factor may be negligibly small or may even be negative. This 
shows the Importance of the bending moment Mp which develops In the 
cylinder as a result of constraining the end. Needless to say, when 
kj becomes negative along any portion of the crack front the problem Is 
no longer a simple crack problem and must be treated as a crack-contact 
problem. In such a case the solution given In this paper Is, of course, 
not valid. The crack contact problem Is highly nonlinear in which the 
additional unknown function L(y) has to be determined by using the Infor- 
mation that the contact Is "smooth" and consequently the total stress 
intensity factor along that portion of the crack front Is zero. Even 
though the problem can be solved by a complicated iterative scheme, 
quite clearly It has no practical value. From a viewpoint of structural 
failure, the important problems In practice are those of crack opening 
not crack closing. 

Finally, some examples showing the distribution of the total stress 
intensity factor kj along the crack front are shown In Fig. 4. Curves 

a and b show kj respectively for an outer and an inner crack for the 

same cylinder and crack dimensions (R/h » 10, L Q /h * 0.4, a/h * l, 
c/a » 1.1). The difference between the two results Is again primarily 
due to bending moment Mp caused by the end constraints. Curve c shows 
the distribution of kj for a longer Internal crack (a/h * 3. c/a ■ 1.1, 

R/h ■ 10, L Q /h - 0.4). The intensification in kj along the crack front 
away from the fixed end observed in curve c may be due to slight cylinder 
bulging as well as to the decrease in the influence of the end constraints. 
Curve d shows kj for an outer crack away from the end (c/a ■ 10, a/h * 1, 

R/h » 10, L Q /h * 0.4) for which the distribution as expected Is nearly 

symmetric. 
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Dimensionless and normalized quantities 


1- 

x » X|/a, y « x 2 /a, z * x^/a, 

„ * C A - d 

1 a 1 1 a 9 

(A. 1 ) 

u » Uj/a, v ■ i,* 2 /a, w * u^/a , 


(A. 2) 

F(x p x 2 ) 

B x - B| i $ y - B 2 , *(x,y) - a5 [ h 

$ 

(A. 3) 

V - °II /E > °yy - ff 22 /E - CT xy - ~T- 

„ » 0, 3 „ „ a 23 

XZ B ’ a yz B » 

(A . 4) 

N . N '> m . N 22 n N I2 

N xx ’ wP N yy * "hT’ N xy TT’ ff lj 

N.. 

- (1-1,2) , 

(A. 5) 

M XX ■ M ll /h2E > \y - M 22 /h2E - \y " M 

12 /h 2 E , 

(A. 6) 

V x ■ V hB > V y * V hB > B * TjfW • 


(A. 7) 

X j ■ 12 (!-v 2 )aVh 2 R 2 , X 4 » 12 (1-v 2 )a 2 /h 2 , k » E/BX 4 . 

(A. 8) 


For dimensions and notation see Figure 1. 


APPENDIX B 


ORIGINAL PAGE ® 
OF POOR QUALITY 


Expressions for stress resultants and rotations. 


N xx (x ' y) 


2 a 4 T..da + 
1 U 


| 2 nj 2 LjjCos3x dB , 


WJ 00 

4 r 8 

N „(x,y) - - 2 a 2 m. 2 T lt da - 2 L n 3 2 cos3x d3, 

yy j i j u i 5 u 


oo 

N xy (x,y) * S a 3 m j .T Jj .da 


00 
f 8 


2 3n.L, .sin 3 x d3, 
5 J *J 


M xx (x ’ y) 


/ 4 

2 (m. 2 -va 2 )T,.da - 
1 J 

-00 


r a 


2 (3 2 -vn?)L-. cos3x d£ 
5 J Z J 


( 1 -v) !a 

IF" 


j . (l-v)a 
ar l T 3 da + "TiP 


3s 


jLjC 0 s 3 x d3 , 


V (x ’ y) 


2 (vm. 2 -a 2 )T,,dce - 
1 J 


r 8 

2 (v3 2 -n. 2 )L,.cos3x d£ 
5 J 

o 


( 1 -v) ia 
" ' h‘P““ 


( 1 -v)a 


or l T 3 d “ • Tu’f 


Sj 3 cos 3 x d3 , 


00 00 
f If f 8 

M xy (x,y) » -I ( 1 -v) 2 am. ? 2 j d“ - ( 1 -v) j 2 3 njL 2 jSin 3 x 


df 


( 1 -v)a 

2 hP 


' <“ 2+r , 2 >V“ + ^ 


(S| 2 +3 2 )LjSinf 


V x (x,y) 


<HA 4 f 4 


2 p.m.T-.da 
J 1 K J J 2 j 


kHA 4 


— j 2 3q.L,.sIn3x d3 
a J 5 J 2 J 


- 1 


0 00 

aTjda + j SjLj sin 3 x d3 , 


(B.l) 
(B.2) . 
(B.3) 

<B.M 

(B.5) 

d3, (B. 6 ) 
(B,7) 
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V y (x,y) 


I«hX 4 


f” 4 


E apjT^jda + 


KhX 4 


00 

r 8 

jj q.n. L- .cosgx dg 
5 J J ^ 
o 


w 

fjTjda - | gL^cosgx dg , 


6 x ( X ,y) - 


00 

f 4 


4 j 8 

E mjTjjdc - Si- j £ BL 2j . singx dB 
0 


-i 


aTjda + 


6 y (x,y) - - 


singx dg , 

oc 

hX 4 


f 4 

E aT„ .da + 

1 2j a 

-00 


f 8 


E n.L 0 . cosgx dg 
5 J Z J 


O 

-I 


r.T-da - 0 B- cosgx dg . 

I 5 Z 


where 


T ij { x»y» a) *271^-7^ 


' ^ j ( a ) iTij x- i ay 


, R. (a) m.x-iay 
T 2j (x - y ’“ ) * 27hXT^7T e 

9 X. 2 n 2 n.y 

L lj (y > 3) ■|T^V^ s j (s)e J ■ 


4j (y ' e) = t hf T wijT-r s j ^ s)e J 


n.y 


LjCy ,3) » -■ (1 7T ~ V - B 2 (g)e Szy 


> 6 . 8 ) 


(B.9) 


(B, 10) 
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Stresses in the uncracked pressurized shell with a fixed end 

in this simple axisymmetric case by assuming that the ends of the 
cylinder is closed and by observing that 


Uj - 0, 3j * 0, e n * u 3 /R, N 22 - pR/2 , 


(C.l) 


following the notation of [18] we can write 

Eh 


N 


22 l-v z dx 


d 2 u. 


Fh du 9 U o 

Eh - V 4 ) . H 


U 3 du 2 

1 1 ~ 1 -v z ^ IT + V dx^ ’ 


m 22 - - D 5^7 - M n * * m 22- d - nnW - 


(C2a,b) 


(C3a-c) 


d M 2 2 j 

djfcT - + ST N n - p 


(C.4) 


where p is the internal pressure and x 2 is the only independent variable. 
From (C.l)-(C.4) it may be shown that 


N . - Eh v£R 
N 11 R u 3 2 > 


(C.5) 


d^u^ 

dx 2 4 


+ A6 4 u 3 = - . (-”<x 2 <0), S 4 - jjpfj- . 


(C.6) 


Solution of (C.6) which is bounded at x 2 =-°° and sat ' s ^' es 


u 3 (0) = 0 , u 3 (0) * 0 


(C.7) 


becomes 
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U 3 ^ X 2 ) “ Ce 2 ( cos 3><2 " s i n 3^2 ) ~ 1 ] , (x 2 <0) . 

From (C.8), (C.5) and (C.3) it then follows that 

2“V ^ x 2 

Ni i (x 2 ) » - ye (cosex 2 -sln3x 2 ) + pR , 

Mj i (x 2 ) 3 -T^ g - 2 — ■ ■ p e (cos3x 2 +sin3x 2 ) , (-«<x 2 <0) . 


(C.8) 

(C.9) 

(C.10) 
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Table 1. Membrane component of the normalized stress Intensity factor k m In 
a pressurized cylinder with a fixed end containing an axial through 
crack, v»0.3. 


— 1 



k n,<- d > 

R/h 

\/a 

a /\ 

1 

1.1 

1.5 

2 

10 

. 

1 

1.1 

1.5 

2 

10 

5 

1 

2 

3 

10 

2.933 

5.75** 

9.16** 

*♦3.555 

0.151 
0.28** 
0. **90 
2.853 

0.3**6 

0.839 

1.500 

*♦.703 

0.572 
1 . 3*40 

2.015 

*♦.912 

1.286 

1.6**5 

2.067 

*♦.923 

0.535 

1.082 

1.583 

*♦.422 

0.575 

1.158 

1.687 

4.605 

0.734 

1.425 

1.985 

4.850 

0.915 

1.613 

2.081 

4.929 

1.284 

1.645 

2.067 

4.953 

10 

1 

2 

3 

10 

J§| 

0.119 

0.181 

0.277 

1.627 

0.235 
0.**77 
0. 8**0 
3. **67 

0.36A 

0.81** 

1.338 

3.605 

I.168 

1.331 

1.589 

3.620 

0.355 

0.713 

1.067 

3.139 

0.360 

0.765 

1.142 

3.314 

0.483 

0.967 

1.402 

3.551 

0.611 

1.167 

1.573 

3.602 

1.158 

1.332 

1.588 

3.620 

25 

1 

2 

3 

10 

1.622 
2. **23 
3-397 
13.151 

0.102 

0.128 

0.165 

0.860 

0.172 

0.275 

0.**2** 

2.218 

0.236" 
0.****8 
0.728 
2. **80 

I.063 

1 . 1**7 
1.259 

2. **69 

0.232 

0.431 

0.650 

2.031 

0.247 

0.463 

0.697 

2.160 

0.303 

0.590 

0.882 

2.429 

0.373 

0.743 

1.074 

2.466 

T.081 

1.142 

1.260 

2.469 


| 

1.345 

1.598 

1.956 

5.563 

0.093 

0.100 

0.111 

0.269 

o. i Ao 
0.169 
0.213 
0.839 

0.167 

0.232 

0.325 

1.336 

0.685 
1.0**5 
1.107 
1 . 5**2 

0.159 

0.229 

0.321 

1.049 

0.166 

0.242 

0.343 

1.126 

0.188“"' 

0.296 

0.434 

1.380 

0.215 

0.367 

0.548 

1.536 

0.744 

1.063 

1.097 

1.542 


| 


0.092 

0.095 

0.101 

0.184 

0. 1 3** 
0.150 
0.173 
0.507 

0.15** 

0.190 

0.2**1 

0.87** 

0.1*96 

0.885 

1.058 

1.308 

0.146 

0.184 

0.238 

0.745 

0.150 

0.193 

0.253 

0.801 

0.165 

0.227 

0.311 

1.007 

O.lTT' 

0.272 

0.388 

1.199 

0.546 

0.937 

1.070 

1.308 
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able 2. Bending component of the normalized stress intensity factor kb In a pressur- 
ized cylinder with a fixed end containing an axial through crack, v * 0.3. 


/h 


k b t- b > 

k b t-d) 

Va 

a /i\ 

1 

1.1 

1.5 

2 

10 

1 

1.1 

1.5 

2 

10 


1 

-20.927 

-0.312 

-0.085 

0.030 

0.115 

0.051 

0.061 

0.100 

0.136 

0.115 


2 

-14.250 

-0.183 

0.069 

0.220 

0.185 

0.151 

0.160 

0,196 

0.211 

0.185 


3 

-11.925 

-0.196 

0.178 

0.266 

0.188 

0.180 

0.189 

0.220 

0.220 

0.189 


10 

-6.839 

-0,284 

-0.582 

-0.869 

-0.935 

-0.491 

-0,593 

-0.808 

-0.920 

-0.933 


1 

-23.036 

-0.316 

-0.134 

-0.051 

0.030 

-0.014 

-0.007 

0.029 

0.069 

0.076 

0 

2 

-13.258 

-0.170 

-0.019 

0.116 

0.147 

0.095 

0.103 

0.138 

0.167 

0.146 


3 

-9.371 

-0.167 

0.077 

0.223 

O.I83 

0.143 

0.153 

0.191 

0.208 

0. 183 


10 

-5.935 

-0.221 

0.025 

-0.144 

-0.199 

-0.029 

-0.046 

-0.111 

-0.181 

-0; 199 


1 

-25. .930 

-0.326 

-0.180 

-0.130 

0.086 

-0.092 

-0.086 

-0.057 

-0.021 

0.073 

5 

2 

-12.693 

-0.165 

-0.091 

0.000 

0.090 

0.023 

0.029 

0.063 

0.098 

0.090 

3 

-8.235 

-0.144 

-0.027 

0.098 

0.138 

0.079 

0.087 

0.123 

0.154 

0.137 


10 

-4.065 

-0.194 

0.256 

0.210 

0.148 

0.159 

0.172 

0.199 

0.174 

0.148 


1 

-28.812 

-0.339 

-0.227 

-0.210 

0.069 

-0.179 

-0.178 

-0. 164 

-0.142 

0.077 

DO 

2 

-12.677 

-0.171 

-0.160 

-0.118 

0.081 

-0.082 

-0.079 

-0.052 

-0.022 

0.070 

3 

-7.162 

-0.140 

-0.116 

-0.052 

0.071 

-0.023 

-0.018 

0.020 

0.048 

0.066 


10 

-1.933 

-0.134 

0.104 

0.243 

0.197 

0.134 

0.147 

0.234 

0.235 

0.198 


1* 

-29.841 

-0.342 

-0.242 

-0.236 

0.019 

-0.209 

-0.209 

-0.202 

-0.187 

0.036 

30 

2 

-12.761 

-0.175 

-0.184 

-0.160 

0.081 

-0.127 

-0.125 

-0.108 

-0.082 

0.082 

3 

-6.911 

-0.144 

-0.148 

-0.106 

0.077 

-0.073 

-0.070 

-0.047 

-0.013 

0.066 

, 

10 

-1.418 

-0.115 

0.018 

0.153 

0.163 

0.091 

0.101 

0.145 

0.180 

0.163 
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Table 3. Normalized stress Intensity factors In a cylindrical shell 
with a fixed end containing an axial semielllptlc surface 
crack and subjected to Internal pressure, R/h - 5. 


Outer Crack, R/h ■ 5, v » 0.3 


a 

£ 

L - 

o 

0.2h 

L ■ 
0 

■ 0.4h 

L o " 

0.6h “ 

L - 0.8h 
0 

h 

a 

i^oT" 

k b (0> 

k t (0) 

■\ToT 

k t (0) 

V oT 

l^ToT 

k b (°) 


1 

.30** 

-.0679 

.185 

-.0367 

.0895 

-.0108 

.0273 

.0009 


1.1 

.331 

-.0490 

.202 

-.0276 

.0984 

-.0087 

.0300 

.0003 

1 

1.5 

.446 

.0160 

.277 

.0064 

.137 

.0009 

.0417 

-.0008 


2 

.585 

.0639 

.368 

.0338 

.183 

.0099 

.0559 

-.0009 


10 

.918 

-.0018 

.588 

-.0011 

.295 

-.0004 

.0898 

.000 


1 

.629 

.0671 

.444 

.0382 

,246 

.0132 

.0818 

.0000 


1.1 

. 680 

.0778 

.484 

.0459 

.270 

.0168 

.0905 

.0007 

2 

1.5 

.856 

.0895 

.621 

.0580 

.355 

.0248 

.121 

.0030 


2 

.979 

.0641 

.718 

.0441 

.415 

.0206 

.143 

.0033 


10 

• 993 

.000 

.737 

.000 

.430 

.000 

.147 

.000 


1 

75to 

.0906 

7673 

.0605 

75To 

.0273 

.150 

.oo4i 


1.1 

.926 

.0859 

.713 

.0590 

.439 

.0279 

.162 

.0048 

3 

1.5 

1.040 

.0489 

.812 

,0373 

.511 

.0 10(> 

.193 

.0051 


2 

1.050 

.0107 

.835 

.0099 

.532 

.0067 

.202 

.0022 


10 

1.020 

.000 

.817 

.000 

.523 

.000 

.148 

.000 


l 

1.060 

-.0020 

.935 

-.0014 

.705 

-.0005 


.0002 


l.l 

1.060 

-.0009 

.935 

-.0008 

.707 

-.0004 

.339 

.000 

10 

1.5 

1.060 

.0001 

.941 

.000 

.714 

.000 

.344 

.000 


2 

1.070 

.000 

.944 

.000 

.720 

.000 

.348 

.000 


10 

1.070 

.000 

.952 

.000 

.732 

.000 

.357 

.000 



* 

Inner Crack, 

R/h - ‘ 

>, v ■ 0.3 




1 

.299 

.0670 

.177 

I .0351 

.0849 

.0099 

.0265 

-.0010 


1 . 1 

• 325 

.0482 

.193 


.0928 

.0080 

.0290 

-.0005 

1 

1.5 

.437 

-.0157 

mim 

-.0059 

.127 

-.0006 

.0399 

.0008 


2 

.572 

-.0625 

.346 

-.0311 

.168 

-.0082 

.0529 

.0014 


10 

.898 

.0018 

.551 


.270 

.0003 

.0844 

.000 


1 

• 6 1 2 

-.0651 

V4l0 

-.417 



.0728 

.0011 


1.1 

.66: 

-.0756 

.445 

-.0411 

.236 

-.0128 

.0797 

.0006 

2 

1 . 5 

.833 

-.0869 

.567 

-.0517 

.305 

-.0190 

.104 

-.0009 


2 

.952 

-.0621 



.355 

-.0157 

.121 

-.0015 


10 

. 968 

.000 

.676 

.000 

.370 

.000 

.126 

.000 


I 

r.857 

-.0876 

7513 


73 W" 

-.0202 

.123 

-.0013 


1.1 

.902 

-.0832 

.649 


.369 

-.0208 

.132 

-.0020 

3 

1 .5 

1.010 

-.0475 

.739 

-.0331 

.427 

-.0156 

.154 

-.0028 


2 

1.030 

-.0105 

.762 

-.0088 

.445 

-.0052 

.161 

-.0014 


10 

1. 000 

.000 

.750 

.000 

.442 

.000 

.160 

.000 


1 

1.050 

-.0018 

.896 

.ooon 

TSzT- 

-.0003 


-.0005 


1 . 1 

1.050 

.0009 

.898 

.0005 

• 6<;2 

.000 

.272 

-.0002 

10 

1 .5 


.000 

.906 


.644 

.000 

.281 

.000 


2 

1.060 


.911 


.652 

.000 

,286 

.000 


10 

1.060 

.000 | 

.920 

.000 1 

.665 

.000 

.294 

.000 
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Table 4. Normalized stress intensity factors Ina cyl indr! cal shell 
with a fixed end containing an axial semielliptic surface 
crack and subjected to Internal pressure, R/h - 10. 

Outer Crack, R/h ■ 10, v - 0.3 


1.1 

.233 

-.116 

.141 

-.0605 

.0676 

-.0168 

.0205 

.0021 

1.5 

.30J4 

-.0535 

.187 

-.0299 

.0907 

-.0093 

.0274 

.0005 

2 

.400 

.0052 

.248 

.0016 

.121 

-.0001 

.0366 

-.0004 


.885 

-.0013 

.559 

-.0006 

.276 

-.0002 

.0829 

.000 

1 

.431 

,0057 

. 300" 

.0006 

.161 ” 

-.0012 

.0515 

-.0008 

1.1 

.470 

,0234 

.330 

.0117 

.178 

.0030 

.0571 

-.0006 

1.5 

.628 

.0721 

.448 

.0447 

.247 

.0174 

.0800 

.0011 

2 

.784 

.0873 

.566 

.0574 

.316 

.0244 

.103 

.0026 

10 

.938 

.0002 

.687 

.0001 

.389 

.000 

.127 

.000 

J 


.0738 


.0472 

. 28 / 

.0195 

“.0982 

.0019 

1.1 

.699 

.0822 

.531 

.0541 

• 315 

.0234 

.109 

.0027 

1.5 

. 866 

.0853 

.670 

.0607 

.408 

.0298 

.144 

.0053 

2 

.972 

.0547 

.760 

.0412 

.468 

.0220 

.167 

.0047 

10 

.969 

.000 

.765 

.000 

.290 

.000 

.170 

.000 

1 

1,030 

-.0003 

.921 

r Too 16 

.705 

.0038 

.337 

.0029 

1.1 

1.030 

-.0032 

.916 

-.0012 

.704 

.0012 

• 339 

.0016 

1.5 

1.010 

-,0013 

.908 

-.0012 

.703 

-.0008 

.340 

-.0002 

2 

1.010 

.0002 

.911 

.000 

.706 

.000 

.343 

.000 

10 

1.020 

.000 

.917 

.000 

.716 

.000 

.350 

.000 


Inner Crack, R/h ■ 10, v ■ 0.3 


71 1.0187 -. 

r*0 ( AAA 1 



t 
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Table 5. Normalized stress intensity factors In a cylindrical 

shell with a fixed end containing an axial semielUptlc 
surface crack and subjected to internal pressure, R/h * 25. 





Outer 

Crack, 

R/h » 25, v * 

0.3 



a 

£ 

L o “ 

0.2h 

L o- 

0.*»h 

L * 
0 

0,6h 

L - 
0 

0.8h 

h 

a 

yoT 

w 

k t (0) 

V°> 

W 

w 

k t (o) 


1 

1 

1.1 

1.5 

2 

10 

.161 

.168 

.20** 

.255 

.853 

-.201 
-.189 
#, i**o 
-.0860 
.0**03 

.0957 

.101 

.12** 

.157 

.53** 

-.102 

-.0972 

-.07**7 

-.0**69 

.0225 

.0456 

.0482 

.0598 

.0759 

.261 

-T 0265 

-.0260 

-.0214 

-.0140 

.0069 

.0138 

.0146 

.0180 

.0228 

.0776 

.0045 

.0040 

.0023 

.0012 

-.0005 

2 

1 

1.1 

1.5 

2 

10 

.27** 

.297 

.397 

.52** 

.908 

-.0909 
-.0721 
-.0053 
. 0**88 
-.0032 

.188 

.205 

.280 

.372 

.656 

-.057** 

-,0**66 

-.0054 

.0306 

-.0022 

.0991 

.109 

.150 

.202 

.362 

-.0219 

-.0185 

-.0031 

.0119 

-.0009 

.0312 

.0342 

.0473 

.0637 

.114 

TTooTT" 

-.0013 

-.0007 

.0006 

.000 

3 

1 

1.1 

1.5 

2 

10 

.*♦11 

.*♦49 

,601 

.757 

.936 

-.0095 

.01*10 

.0669 

.0879 

.0002 

.306 

.336 

.*♦58 

.582 

.729 

-.0057 

.0069 

.0454 

.0628 

.0001 

.174 

.193 

.268 

.345 

.439 

-.0039 

.0015 

.0202 

.0303 

.000 

.0568 

.0631 

.0887 

.116 

.148 

-.0011 

-.0005 

.0024 

,0047 

.000 

10 

1 

1.1 

1.5 

2 

i_UL_ 

1.000 

1.020 

1.010 

.991 

.995 


. 89 ** 

.908 

.915 

.898 

.899 

.0420 
.0317 
.0031 
-.0028 
. .000 

.675 
.691 
.708 
.698 
_^701 

.0299 

.0239 

.0043 

-.0017 

.000 

.306 
.317 
.332 
.329 
• 3”?2 

.0108 

.0094 

.0028 

-.0004 

.000 




Inner 

Crack, 

R/h * 25, V - 0.3 



1 

1 

.160 

.168 

.203 

.253 

.8**7 

.200 

.189 

.l**0 

.0855 

-.0**00 

.09**8 

.0998 

.122 

.15** 

.523 

.101 

.0960 

.0734 

.0459 

-.0218 

.0451 

.0476 

.0588 

.0743 

.253 

.0257 

.0252 

.0206 

.0133 

-.0065 

.0138 

.0145 

.0178 

.0225 

BBQi 

2 

1 

1.1 

1.5 

2 

10 

.272 

.29** 

.39** 

.518 

.898 

.090** 
.0716 
.0053 
- . 0**83 
.0032 

.185 

.201 

.272 

.360 

.632 

HRK ; : jj|| 

t ; j| 
1 & 

.0960 

.105 

.143 

.191 

.339 

.0211 

.0176 

.0030 

-.0107 

.0008 

11 

.0010 

.0011 

.0007 

-.0003 

.000 

3 

1 

1.1 

1.5 

2 

10 

.**07 
.444 
.59** 
• 7**7 
.925 

.00**6 

-.0137 

-.0661 

-.0867 

-.0002 

.296 

.32** 

.*♦38 

.555 

.696 

.0060 

-.0062 

-.0431 

-.0593 

-.0001 

.164 

.180 

.247 

.316 

.401 

.0042 

-.0009 

-.0177 

-.0263 

.000 

i 

.0013 

-.0008 

-.0015 

-.0031 

.000 

10 

1 

1.1 

1.5 

2 

10 

.991 

1.000 

1.000 

.982 

.981 

-.03**! 
-.0017 
.003** 
— ,.QQ.Q 

\#*r 
.86** 
.87** 
.860 
.Ml . , 

-.0296 

-.0032 

.0026 

■iQ.QQ., 

“.592 

.606 

.625 

.620 

.626 

-“.0246 

-.0199 

-.0040 

.0013 

.239 

.248 
.261 
.261 
• 266 

-.0069 
-.0062 
-.0021 
.0002 
>.QQ.O. . , 
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Table 7. Normalized stress Intensity factors In a cylindrical 

shell with a fixed end containing an axial semlelllptlc 
surface crack and subjected to Internal pressure, R/h ■ 200. 





Outer 

Crack, 

R/h » 200, V * 

0.3 



a 

c 

L * 

0 

Q.2h 

L - 

o 

0.4h 

L o- 

Oh 

L o- 

Oh 

h 

a 




J£To) 

k t (or 

yo7— 

k t ToT 


1 

1 

1.1 

1.5 

2 

10 

.124 

.126 

.133 

.142 

.426 

-.297 

-.293 

-.273 

-.247 

.0261 

.0734 

.0749 

.0802 

.0867 

.264 

-.150 

-.149 

-.144 

-.132 

.0143 

.0349 

.0357 

.0385 

,0418 

.128 

-.0385 

-.0395 

-.0401 

-.0380 

.0043 

.0106 

,0108 

.0116 

.0126 

.0380 

.0069 

.0064 

-.0227 

.0041 

-.0004 

2 

1 

1.1 

1.5 

2 

10 

.152 

.156 

.177 

.209 

.797 

-. 26 } 

-.253 

-.212 

-.162 

.0801 

.103 

.107 

.124 

.147 

.568 

■^.T 62 

-.158 

-.136 

-.106 

.0532 

.0535 

.0558 

.0654 

.0783 

.306 

-.0595 

-.0591 

-.0533 

-.0426 

.0221 

.0168 

.0175 

.0204 

.0243 

.0948 

-.0018 

-.0023 

-.0031 

-.0028 

.0017 

3 

1 

1.1 

1.5 

2 

10 

.183 

.192 

.234 

.295 

.942 

-.216 

-.202 

-.144 

-.0821 

.0364 

.134 

.142 

.176 

.223 

.721 

- .1 48 
-.139 
-.103 
-.0596 
.0266 

.0747“ 

.0794 

.0998 

.128 

.420 

-.0637 

-.06l4 

-.0474 

-.0284 

.0129 

.0239 

.0254 

.0320 

.0410 

.135 

-.0065 

-.0067 

-.0060 

-.0039 

.0019 

10 

1 

11 

1.5 

2 

10 

.504 

.549 

.719 

.873 

.972 

-.0314 

.0482 

.0876 

.0880 

.000 

.444 

.485 

.639 

.780 

.874 

.0242 

.0387 

.0744 

.0767 

.oca 

"TITS 

.349 

.468 

.577 

.652 

.0131 

.0227 

.0486 

.0525 

.000 

.124“ 

,137 

.190 

.238 

.273 

.0026 

.0054 

.0142 

.0168 

.000 




Inner 

Crack, 

R/h * 200, v * 

0.3 



1 

"T~ 

1.1 

1.5 

2 

10 

. 124 H 
.126 
.132 
.142 
.425 

.297 

.293 

.273 

.247 

-.0261 

.0734 

.0748 

.0801 

.0866 

.263 

.149 

.149 

.143 

.132 

-.0142 

.0349 

.0357 

.0385 

.0417 

.127 

.0364 

.0394 

.0400 

.0378 

-.0042 

.0106 

.0108 

.0116 

.0126 

.0378 

-.0069 

-.0064 

.0051 

-.0042 

.0004 

2 

1 

1.1 

1.5 

2 

10 

.151 

.156 

.177 

.209 

.795 

.263 

.253 

.211 

.162 

-.0799 

.103 

.107 

.123 

.146 

.563 

.162 " 
.157 
.135 
.105 
-.0527 

.0533 

.0556 

.0650 

.0777 

.302 

.0591 

.0587 

.0527 

.0420 

-.0215 

.0167 

.0174 

.0203 

.0242 

.0934 

.0017 

.0022 

.0029 

.0027 

-.0015 

3 

1 

1.1 

1.5 

2 

10 

.183 

.191 

.233 

.294 

.939 

.216 

.202 

.144 

.0819 

-.0363 

.134 

.141 

.174 

.221 

.713 

.147 

.139 

.102 

.0589 

-.0262 

.0741 

.0786 

.0985 

.126 

.410 

. 0629 ” 

.0605 

.0465 

.0276 

-.0124 

.0237 

.0252 

.0315 

.0403 

.132 

.OO 63 

.0065 

.0057 

.0036 

-.0017 

10 

1 

1.1 

1.5 

2 

10 

.502 

.546 

.715 

.868 

.967 

.0312 

-.0479 

-.0871 

-.0876 

.000 

.436 

.475 

.625 

.762 

.855 

-.0233 

-.0376 

-.0725 

-.0746 

.000 

• 303 
.332 
.441 
.542 
.614 

-.0118 

-.0210 

-.0451 

-.0484 

.000 

.114 

.126 

.170 

.212 

.243 

-.0018 

-.0043 

-.0118 

-.0137 

.000 


Table 8. Distribution of the normalized stress Intensity factors 
along the crack front In a cylindrical shell containing 
an axial seml-el 1 1 ptfc surface crack and subjected to 
Internal pressure, x ■ (x«+c)/a, R/h ■ 10, a/h ■ 1, c/a ■ 
1 . 1 . 



0.2 

0.4 

0.6 

“ 0.8 ~~ 


m 

m 

1 m 

mm 

n 

mm 

ra 

k b 

X 

Outer Crack, R/h 

■ 10 , a/h ■ 1 , 

c/a ■ 1 

1 


mmm 

WiKfcW 


.0615 

PW 

.0310 

jp ;■ 


-.0212 

.828 

.106 

-.247 


-.145 

.0339 

Ojire 

WSSm 

-.0128 

.688 

.127 

-.231 

.0806 


.0394 

-.0442 

Ell 

-.0055 

.516 

.153 

-.203 


-.102 

.0467 

-.0318 

Mm 

-.0005 

.319 

.183 

-.169 

.113 

-.0837 

.0556 

-.0229 

.0168 

.0019 

.108 

.216 

-.133 

.131 

-.0676 

.0635 

-.0181 

.0193 

.0026 

0 

.233 

-.116 

.141 

-.0605 

.0676 

-.0168 

.0205 

, 'frill 

-.108 

.248 

-.0990 

.149 

-.0539 

.0716 

-.0161 

.0216 


-.319 

.275 

-.0687 

.164 

-.0420 

.0785 

-.0153 

.0236 


-.516 

.293 

-.0433 

.175 

-.0316 

.0829 

-.0144 

.0257 

158® 

-.688 

.298 

-.0231 

.179 

-.0224 

.0846 

-.0126 

.0267 

-.0034 

-.828 

.288 

-.0085 

.176 

-.0147 

.0840 

-.0103 

.0259 


-.929 

.258 


.167 


.0818 

-.0082 

.0251 

-.0031 


Inner Crack, R/h 

- 10, a/h * 1 , 

c/a ■ 1 

1 


.929 


.247 

.0619 

.169 

.0310 

.0794 

.0094 

.0210 

.828 

.106 

.246 

.0683 

.144 

.0334 

.0582 

.0104 

.0125 

.688 

.127 

.230 

.0796 

.121 

.0386 


.0124 

.0052 

.516 

.152 

.202 

.0941 

.100 

.0454 

.0306 

.0144 

.0003 

.319 

.182 

.167 

.110 

.0818 


.0217 

.0164 

-.0022 

.108 

.2)4 

.132 

.128 

.0658 

.0615 

.0170 

.0189 

-.0028 

0 

.230 

.114 

.137 

.0587 

.0655 

.0158 

.0201 

-.0023 

-.108 

.246 

.0980 

.145 

.0522 

.0694 

.0151 

.0212 

-.0013 

-.319 

.273 

.0680 

.160 

.0407 

.0761 

.0144 

.0231 

.0008 

-.516 

.291 

.0428 

.171 

.0306 

.0805 

.0136 

.0251 

.0024 

-.688 

.297 


.175 

.0217 

.0821 

.0119 

.0260 

.0032 

-.828 

.286 

.0084 

.173 

.0142 

.0817 

.0097 

.0251 

.0032 

-.929 

.257 

.0004 

.165 

.0085 

.0796 

.0077 

.0242 

.0029 


PAGE te 

OF POOR QUALITY 


-4o- 











ORIGINAL PAGE !9| 

OF POOR QUALITY 

Table 9. Distribution of the normalized stress Intensity factors 
along the crack front in a cylindrical shell containing 
an axial semi-elliptic surface crack and subjected to 
internal pressure, x - (xj+cJ/a, R/h * 10, a/h * I, 


L o' h 

0.2 

0.4 

0.6 

0 

.8 


k t 

k b 

k t 

k b 

k t 

k 

K b 

k t 

k b 

X 

Outer Crack, R/h « 10, a/h * 1 

, c/a - 
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IT 

. Tol 

- .0858 

.124 

- .0555 

.0588 

-'.'6452 

.522T“ 


.828 

.202 

-.0813 

.143 

-.0456 

.0755 

-.0180 

.0246 

-.0038 

.688 

.245 

- .0668 

.165 

-.0338 

.0845 

-.0116 

.0277 

-.0012 

.516 

.291 

-.0469 

.190 

-,0222 

.0951 

-.0062 

.0302 

.0002 

.319 

.337 

-.0252 

.214 

-.0117 

.106 

-.0028 

.0324 

.0005 

.108 

.381 

-.0044 

.238 

-.0025 

.117 

-,0008 

.0352 

.000 

0 

.400 

.0052 

.248 

.0016 

.121 

-.0001 

. 0366 

-.0004 

-.108 

.418 

.0141 

.258 

.0055 

.126 

.0006 

.0378 

-.0007 

-.319 

.443 

.0292 

,272 

.0126 

.132 

.0025 

.0399 

-.0008 

-.516 

.454 

.0404 

.278 

.0190 

.135 

.0052 

.0421 

-.0003 

- . 688 

.446 

.0476 

.277 

.0247 

.134 

.0085 

.0429 

.0009 

- . 828 

.419 

.0503 

.266 

.0294 

.131 

.0119 

.0411 

.0026 

--929 

.369 

.0480 

.249 

.0327 

.126 

,0153 

.0396 

.0041 


Outer Crack, R/h 

* 10, a/h ■ 1 , 

c/a * 10 


~ W 

.592 

.0022 

.429 

.0013 

.230 

.0005 

.0742 

.0001 

.828 

.693 

.0020 

.468 

.0010 

.241 

.0003 

.0774 

.000 

.688 

.770 

.0013 

.503 

.0006 

.253 

.0001 

,0821 

.000 

.516 

.827 

.0006 

.531 

.0002 

.265 

.000 

.0832 

.000 

.319 

.865 

-.0002 

.549 

-.0002 

.273 

-.0001 

.0825 

.000 

.108 

.883 

-.0009 

.558 

-.0005 

.276 

-.0002 

.0828 

.000 

0 

.885 

-.0013 

.559 

-.0006 

.276 

-.0002 

.0829 

.000 

-.108 

.881 

-.0016 

.558 

-.0008 

.276 

-.0002 

.0827 

.000 

-.319 

.860 

-.0020 

.547 

-.0010 

.272 

-.0003 

.0822 

.000 

-.516 

.820 

-.0024 

.527 

-.0013 

.263 

-.0004 

.0827 

.000 

-.688 

.762 

-.0026 

.498 

-.0015 

.251 

-.0006 

.0815 

.000 

-.828 

.684 

-.0024 

.463 

-.0016 

.239 

-.0007 

.0768 

-.0002 

-.929 

.583 

-.0013 

.424 

-.0017 

.227 

-.0009 

.0736 

-.0002 
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Table 10. Distribution of the normalized stress Intensity factors 
along the crack front in a cylindrical shell containing 
an axial semi -el 1 iptlc surface crack and subjected to 
internal pressure, x ■ (xj+c)/a, R/h »• 10 , c/a * 1.1. 
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b 
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X 

Inner Crack, R/h 
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.T0h 
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.0571 

H563T" 

'.0236 

.0l$l 

.828 

.159 

.143 

.126 

.0912 

.0811 

.0414 

.0317 

.0106 

.688 

.244 

.0829 

.191 

.0476 

.117 

.0187 

.0439 

.0031 

.516 

.362 

.0169 

.275 

.0069 

.162 

.0007 

.0581 

-.0006 

.319 

.498 

-.0373 

.369 

-.0242 

.212 

-.0104 

.0729 

-.0014 

.108 

.629 

-.0711 

.458 

-.0436 

.258 

-.0166 

.0879 

-.0010 

0 

.685 

-.0804 

.495 

-.0494 

.277 

-.0187 

.0944 

-.0010 

-.108 

■ 731 

- . 0850 

.526 

-.0529 

.293 

-.0205 

,0993 

-.0014 

-.319 

.787 

-.0840 

.563 

-.0547 

.313 

-.0233 

.106 

-.0033 

-.516 

.794 

-.0745 

.566 

-.0517 

.316 

-.0250 

.109 

-.0053 

-.688 

.755 

-.0619 

.538 

- . 0466 

.302 

-.0254 

.107 

-.0070 

-.828 

.676 

-.0488 

.483 

-.0401 

.277 

-.0241 

.0989 

-.0080 

-.929 

.559 

-.0364 

.409 

-.0326 

.243 

-.0217 

.0894 

-.0081 


Inner Crack, R/h * 10, a/h » 10, 
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.1 


• 929 


.0164 

.176 

.0091 

.123 

.0023 

.0598 

-.0005 

.828 

.431 

-.0519 

.334 

-.0411 

.226 

-.0254 

.103 

-.0094 

.688 

.688 

-.0791 

.544 

-.0584 

.366 

-.0332 

.160 

-.0105 

.516 

.894 

-.0547 

.726 

-.0410 

.494 

-.0235 

.212 

-.0071 

.319 

.996 

-.0187 

.835 

-.0163 

.579 

-.0110 

.244 

-.0039 

.108 

1.020 

.000 

.874 

-.0024 

.616 

-.0037 

.260 

-.0021 

0 

1.010 

.0029 

.874 

.0003 

.619 

-.0018 

.263 

-.0015 

-.108 

1.010 

.0033 

.867 

.0012 

.615 

-.0008 

.261 

-.0010 

-.319 

.965 

.0018 

.822 

.0011 

.581 

.000 

.249 

-.0005 

-.516 

.898 

.0006 

.744 

.0005 

.518 

.000 

.227 

-.0003 

-.688 

.805 

-.0003 

.643 

-.0006 

.438 

-.0006 

.195 

-.0004 

-.828 

.688 

.0001 

.526 

.0005 

.350 

.0004 

.157 

.000 

-.929 

.549 

.0001 

.407 

.0004 

.266 

.0004 

.120 

,000 
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Fig.' 3 Normalized total stress Intensity factor at the maximum pene- 
tration point (L^Lq, X 2 as- c) of a semi -el 1 iptic outer crack in a 
pressurized cylinder with a fixed end, v * 0.3, L 0 /h = 0.4, 




Distribution of the normalized total stress intensity factor along the border 
of a part-through semi -el 1 iptic surface crack in a pressurized cylinder with 
a fixed end, v - 0.3» L 0 /h - OA, P t /h = 10. Curve a: outer crack, a/h = 1, 
c/a = 1.1; curve b: inner crack, a/h = 1, c/a = 1.1; curve c: Inner crack, 
a/h = 3, c/a = 1.1; curve d: outer crack, a/h = 1, c/a = 10. 


